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The continuity equations expressing conservation of string defect charge can be used to find an 
explicit expression for the string velocity field in terms of the order parameter in the case of an 0(n) 
symmetric time-dependent Ginzburg-Landau model. This expression for the velocity is used to find 
the string velocity probability distribution in the case of phase-ordering kinetics for a nonconserved 
order parameter. For long times t after the quench, velocities scale as t -1 / 2 . There is a large velocity 
tail in the distribution corresponding to annihilation of defects which goes as y-( 2d + 2 - n ) f or both 
point and string defects in d spatial dimensions. 

PACS numbers: 05.70.Ln, 64.60.Cn, 64.75.+g, 98.80.Cq 



I. INTRODUCTION 

In recent workEI we discussed how one could use conservation of topological charge to study the statistics of velocities 
of point defects in phase-ordering systems. We were able to identify the appropriate point-defect velocity field in terms 
of the order- parameter field in the context of a <i-dimensional 0(n) symmetric time-dependent Ginzburg-Landau 
(TDGL) model. Using this expression for the velocity field, for point particles (n = d), the probability distribution 
for defect velocities was determined in the case of the late-state phase ordering using the lowest-order approximation 
in the perturbation expansion method developed in Ref.( |^). This analysis is extended here to the case of strings 
defects where n = d — 1. The velocity probability distribution is worked out explicitly for n = d — 1 = 1 and 2 and 
for the wall case n = d — 2 = 1 . These results and the results for all n — d can all be written in the form 

/ 1 \ d/2 mra+D fv^-*( v (d+2)/2 

P( V )= -L UJ >T (tf) i + v 2 /« 2 • (i) 

where the velocities scale with a factor v w L(t)~ x where L(t) « t x l 2 is the characteristic scaling length which grows 
with time t after the quench. The result for P(V) indicates that the probability of finding a defect with a large 
velocity decreases with time. There is a high- velocity tail y-( 2d + 2 - n ) which corresponds to the annihilation of defects 
and defect loops. It is encouraging that this result agrees with the large V tails obtained by BrayS using rather general 
scaling ideas for all n < d. 

The result obtained here for P[V] seems very simple. Does it correspond to experimental observation or the results 
of numerical simulations? Thus far there have been not direct tests. It would seem worthwhile to check the range of 
validity of the defect velocity probability distribution given by Eq. ([!]). 

II. PROBLEM SET UP 

We study an n-component non-conserved order-parameter field ip a (R, t) in (i-spatial dimensions which satisfies the 
TDGL equation 

^ (R ' t) -- r ^S) + ^ (R ' t) (2) 

where F is an effective free-energy functional and T is a constant kinetic coefficient. We assume F is of the 0(n) 
symmetric square-gradient form 
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where c > and V(ip) is chosen to be a degenerate double-well or wine-bottle potential. This model is to be 
supplemented by random, uncorrelated, initial conditions. We assume that there is a rapid temperature quench from 
a high temperature to zero temperature where the noise r\ a in (0) can be set to zero. In the scalar case (n=l) such 
systems order through the growth of domains separated by sharp walls. As time evolves these domains coarsen and 
order grows to progressively longer length scales. In the case of systems with continuous symmetry (n > 1) the 
disordering elemcntsu'EI will depend on n and spatial dimensionality d. Thus, for example, for n — d one has point 
defects (vortices or monopoles) while for n = d — 1 one has vortex lines or string-like objects. For n > d there are no 
stable singular topological objects. ^ p 

The main physics in phase-ordering systemsBia is the interplay between two characteristic lengths, a characteristic 
domain size L(t), which grows with time, and a defect dimension £ (interfacial width, vortex core size, etc.). However, 
at long enough times the single length L(t) dominates, L(t) ^> £, the morphological structure looks self-similar under 
the rescaling of space and time and the order parameter correlation function satisfies the scaling equation 

CV(R, t) = (V7(R, t) ■ $0, t)) = tl$F(x) (4) 

where x = R/L(t), and is the magnitude of ip in the ordered state. The structure factor, the fourier transform of 
C^(R, t), satisfies CL,(q,t) = L d ipQF(Q), where Q = qL is a scaled wavenumber. For pure systems with short-range 
interactions and a nonconserved order parameter the growth law is given by the Lifshitzj-.Cahn- Allen result L « t 1 / 2 
for all n. The two-dimensional XY modeHTJ and the one-dimensional scalar modelolij appear to be interestiaa 
exceptions. For large Q and n < d, due to defects, the structure factor obeys the generalized Porod's Law,E3~E2l 
HQ) ~ Q~ (n+d) - This reflects increasingly weaker singularities in F(x) for small a; as a function of n. In the opposite 
limit, it appears that the large ff behavior can, with proper definition of x, be put in the form F(x) ~ x~ v e~^ x 
where v is a subdominant indexB. The quantities-.discussed above are evaluated at equal times after the quench. In 
the two-time case one again has a scaling lawEiLJ and the on-site correlation function has the form 

< ^(R, t + 1) ■ $R, t) >« L{t)- x , (5) 



for r>i, where A is a nontrivial exponent which has been determined numerically and theoreticallya for a number 
of systems. 

III. DEFECT DYNAMICS FOR POINT PARTICLES 

Since a great deal . is known about order-parameter correlations in phase-ordering systems, attention has turned 
toward the studvE-aEj of the statistics and dynamics of the annihilating defects themselves. The basic idea is that the 
positions of defects are located by the zeros of the order parameter field tp, therefore the charged or signed density 
for point defects is given by 

p(R, t) = <5(Vj(R, t))V(R, t) (6) 
where T>, associated with the change of variables from the set of vortex positions to the field ip, is defined by: 

T>= ^W2...a^^2...^ v pi?/Vl V a'2' ! /v 2 ■ ■ • V M „Vv„ (7) 

where e^ lM2 ... Mjl is the n-dimensional fully anti-symmetric tensor and summation over repeated indices is implied here 
and below. The unsigned defect density, n(R, t), is given by 
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(R,t)=^(R,t))|2?(R,t)| . (8) 



It was shown in Ref.( [!]) that the vortex charge density p for point defects, defined by Eq.(^), satisfies a continuity 
equation of the form 

p = -v M1 KJ (9) 

where the defect velocity field w Ml is given by 
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where V is defined by Eq. (|7|) . Thus one has an explicit expression for the defect velocity field which can be expressed 
strictly in terms of the order parameter and its spatial derivatives. Remember that the TDGL equation of motion 
can be used to express tp Ul in terms of ip and it spatial derivatives. 

The expression given by Eq.(|l0|) for the velocity is very useful because it avoids the problem of having to specify 
the positions of the defects explicitly. The positions are implicitly determined by the zeros of the order-parameter 
field. The practical usefulness of Eq.(|l0|) can be seen by asking the question: In the scaling regime of a phase-ordering 
system with point defects, what is the probability of finding a defect with a velocity V? This probability distribution 
function is defined by 

(n)P(V) = (nS(V-y(iJ))) , (11) 

where V is a reference velocity and n is the unsigned defect density defined by Eq.(||). The calculated P(V) is given 
by Eq.(||) with n = d, _ 

Going further along these lineso we considered the two vortex velocity probability distribution, P[V 1; V 2 , R] which 
gives one the probability of finding the velocity of one defect in the fixed presence of another defect a known distance 
away with a known velocity. Clearly P[Vi, V2, R] contains a tremendous amount of information about the dynamics 
of point defects. The physical results from the calculation of this quantity, carried out in detail for n = d = 2 in REf.( 
p6| , are relatively simple to state. The probability distribution is a function only of the scaled velocities Uj = Vi/v for 
i = 1 or 2, and the scaled separation x = R/L(t). The characteristic velocity v is the same quantity that appears in 
P(V). For a given scaled separation x , the most probable configuration corresponds, as expected, to a state with zero 
total velocity and a nonzero relative velocity only along the axis connecting the vortices: Vi = — V 2 = vx. Moreover 
there is a definite most probable nonzero value for v — v max for a given value of x. The most striking feature of these 
results is that for small x the most probable velocity goes as v max = k/R where R is the unsealed separation between 
the vortices and n = 2.19 in dimensionless units defined in Ref.( ^(f ). The result giving v max inversely proportional 
to R is consistent with overdamped dynamics where the relative velocity of the two vortices is proportional to the 
force which in .turn is the derivative of a potential which is logarithmic in the separation distance. Since there is 
low probabilityE-ffiH of finding like-signed vortices at short distances, our results giving the velocity as a function of 
separation distance should be interpreted in terms of annihilating vortex-antivortex pairs and is in agreement with 
the general scaling ideas proposed by BrayQ. 



IV. CONTINUITY EQUATIONS 



Let us investigate the existence of a local statement of topological charge in the general case of n < d. The first 
step is to introduce the appropriate density for topological charge. Jii the case of point particles the conserved density 
is the charge density given by Eq.(^|). The next obvious extensionEa is to string defects where n = d — 1 and the the 
defect line density is given by 

Psi = 8$)V S1 (12) 

where 

and the s and [i range from 1 to d and the v range from 1 to n and there is summation over repeated indices. The 
generalization of Eqs.([|) and (|l2| ) to all n < d is given by 

p slS2 ... Sd _„ = 8{^)V slS2 ,„ Sd _ n (14) 

where 

It is then straightforward to show that T> slS2 _ Sd _ n itself satisfies a continuity equation given by 

2^SiS 2 ...Sd_„ — Vjtj /siS 2 ...S d _„/ii j (16) 

with the current defined by 
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and 

To obtain the continuity equation for p Sl s 2 ...s d „ rl , we need a second identity. Consider the quantity 

JsiS2...S d - n HlV Ultyu = e siS 2 ...Sd_„MlM2---Mn$M2A I 3---A I r 1 ^Ml^ I/ 
e SiS2...S £ j_„^Ii^2---Mn e i / l' y 2---^r 1 ^Ml 1 /' iy ^M2 VV 2 ■ ■ ■ ^ ^n^P^n ' (l^) 

The key observation is that the right-hand-side of Eq.([l9|) has a factor which can be written in the form: 

£ SlS 2 ...S d -n^l^2 ■••A'n ^/V^M2 ^Va ' ' " ^(ln^"n = E W2...I'„QsiS2...S J „„ (20) 



where Q SlS2 ...s i _ n is determined by multiplying Eq.(f20|) by ej/,^...j/ n and summing over all the i/s. We easily obtain, 
remembering Eq.(p"5|), 

Qsis 2 ...s d „„ — ^sis 2 ...s d _„ • (21) 

Putting this result back into Eq.(M) gives 



1 

{n~ 1)!' 



■7siS2...S^_ n /il V^j^^y — . . V , fi e fl^ 2 .--^„ e ^^ 2 -..^„^'s 1 S 2 ...S d _ 



- — ;ip Vl S Vl , v (n- l)W SlS2 ... Sd _ n 



(n-l)l 



Taking the time derivative of Psis 2 ...s d _„ gives, using Eqs. (10) and (E2J), 



dip v 

and finally we obtain the desired continuity equation 

Ps 1S2 ...s d -„ = V m ^W)J SlS2 ... Sd _„ w ) • (23) 
For the simplest case of point defects (n = d), Eq.(]23]) can be put into the conventional form given by Eq.([|) with 

Jin = -v„V > ( 24 ) 

and the velocity given by Eq.(|l0|). 

Let us turn next to the case of strings where the line density is a vector p Ml and the current is a two-component 
tensor, 

Clearly J Sl/Jl is anti-symmetric in its subscripts. Since we expect the instantaneous velocity to be orthogonal to the 
local orientation of the string, we can define the velocity via 

Ja/3 = V a T> - VpT> a . (26) 
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Dotting the vector T> into this expression gives the result 



Va = -^-JapVp , (27) 



where we have taken advantage of the fact that v and T> are orthogonal: 



■D-v = V a: —J a0 Vp = O . (28) 
D 1 

The velocity field for strings can be written in the form 
for general n. 

Let us check this result and its sign for the simplest case of n = 1, d = 2. The vector V in this case takes the simple 
form 

£Vl = W2 V M2'*/' > ( 30 ) 

9 = i> , (31) 

and 

Wsi = (V^ W2VM2 ^ esiiUl ^ = _ (v|)2^ ' (32) 

Consider a circular loop of string of radius R(t) where the order parameter near the interface formed by the loop is 
given in polar coordinates in the form 

tp(r) = A(r ~ R{t)) . (33) 
where A is an overall constant amplitude. We then need the derivatives 

V^ = Af^ (34) 

and 

= -R(t)A (35) 

to obtain the velocity 

v Sl = -^R(t)r Sl A 2 = R(t)r Sl (36) 

as expected. Typically we will substitute for i\) using the equation of motion and use the defect locating <5-function to 
set 

i> = IW 2 V> (37) 

and 



Using the same ansatz given by Eq.(|33J) and remembering that the expression for the velocity is multiplied by a 
(5-function setting r = R(t), leads to the result: 

^ = - = W) ' (39) 
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and we obtain the Lifshitz-Cahn-AllenEl □ result, 



v ^ = - Tc wf si ' (40) 



which tells us that the circular domain is shrinking and R(t) ~ t 1 / 2 . 
For n = 2 and d = 3 we have explicitly the results reported in Rcf. ( 



V 
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2<*im (fip ui x Vt/v 2 ) (41) 
9 '= e^V^V^ (42) 



and 



u= x g . (43) 



V. BRIEF DISCUSSION OF WALLS 

Let us briefly discuss how things develop as one attempts to go further and increase d — n to 2 and the case of 
walls. We must in this case deal with the quantities 

1~^S1S2 ~\ £ SlS2Ail/J2---/Jn ^V\V 2 ---^n VVl ^Hn^Vn (44) 

and 

JsiSzfll — e SiS2AIlP2---Mn9A 1 2A 1 3---Mn " (4^) 

If we recall the definitions of the velocity field for point defects 

= -v^V (46) 

and for string defects 

Jap = v a T>p - vpV a , (47) 

then it is natural to write for walls 

This expression builds in the antisymmetry of J Sl s 2 m- 

If one restricts the discussion to the physically most relevant case of n — 1, d = 3, then one can use the result 

V SlS2 = f SlS2W V w i/) (49) 

to show that the velocity field is given by 

1 

V 2 ' 

where 



= -^2^182^1*2^3^1^ (50) 



V 2 = (Vipj . (51) 
Further straight-forward manipulation leads to the final result 

- ~ V ^ j, • (52) 



The key result here is that v is orthogonal to X> SljS2 : 



v^V^ss = rjV»iV> tn-L,'*,^ n 3 ip = . (53) 



It remains to be seen if Eq.(|48[) serves as a useful definition of the velocity field for walls with n > 1. Notice that the 
results for a scalar order parameter (n = 1) can all be written in the form of Eq.([52"|) for d = 1,2 and 3. 
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VI. EVALUATION OF P[V] FOR STRINGS 



Our interest here is in ^determining the defect velocity probability distribution, P[V], for strings. Again we use 
the auxiliary field methodcjcj which has been successful in determining the scaling function for the order-parameter 
correlation function in a perturbation theory expansion. We evaluate P here to lowest order in this expansion where 
the auxiliary field can be treated as a gaussian field. The first step in this theory is to express the order parameter 
in terms of an auxiliary field rh. For our purposes here the important result is that near a charge one vortex core the 
order parameter is linear in the auxiliary field 

ip(m) = Am + 0(m 3 ) . (54) 



It is then easy to show that one can replace ip by m in the expression for v given by Eq. ( |29| ) and in the expression for 
the string-charge density 

Patty) = Pa(m) . (55) 
We then want to determine the string-velocity probability distribution 

(\p\)P(V) = {\p(j)\6(V-v($))) , 

= (\p(m)\S(V - v(m))) . (56) 
One can determine P(V) by first evaluating the more general probability distribution 

G(£, b) = (S(m)5(t; - V^Sib- V 2 m)) (57) 



since 



/d n 
d n bH ]Jd%\V(mS(V-t(b,Z))G(£,b) (58) 
, i .. — i 



fl — l u — l 



where 



with 



and 



Tc 1 



^Sl(^) - < e eifJ,lfJ,2---fJ.n e i'lV2---^n^fl 1 1 ifJ 2 2 ■ ■ ■ ^flZ V*0) 



n = (\p(m\) . (61) 

We have assumed that the quench is to zero temperature so that the noise can be set to zero and we can use Eq.([37|). 
The gaussian average determining G(£, b) is relatively straightforward to evaluated and is given by 



1 e 2S 4 1 
G( ^' 6) = (27rS* )«/2 ( 2 ^ 4 )"/ 2 (2nS 2 ) nd / 2eXP 



■ 2S Y,(0 

{J,, I* 



(62) 



where 

1 



5 = -{m z ) S3 V (63) 



77 



§ 2 = ^-{(Vrh) 2 )^L° (64) 
an 
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and 



5 4 = -<(V 2 m) 2 ) - ' 

n S 



(65) 



The quantities Sq, §2, S4 arc determined from the theory for the order parameter correlation function and discussed 
further below. 

The problem then reduces to evaluating the b and £ integrations in the integral given by Eq.([38|) using the result 
for G(£, b) given by Eq.(p2). We proceed by first doing the integration over b. This is facilitated by first defining the 



matrix MY. via 



and 



«„ = YcMY\b v 



(66) 



MY 



(»-l)!23 2 (0 



(67) 



Clearly the quantity T> 2 is important in the development and is discussed in Appendix A. The matrix M" is discussed 
in Appendix B. Then we use the integral representation for the 5-function and find that we can evaluate the integral 
over b in terms of standard displaced gaussian integrals with the results: 



"° p[vl =^/nn«(^ 



exp 



where 



and 



/( = 1 V — 1 



m 



2S& 



f t . V 



d d k 

(2Tr) d ' 



Ak-V 



MO 



MO 



d% 

(27^4 ) n/2 



(68) 



(69) 



(2 7 rS , 4 )"/ 2 



exp 



-S^TcfkakfiMlMl 



(70) 



Sums over a and (3 range from 1 to d, those over v's from 1 to n. We then have the apparently Gaussian integral over 
k: 



J(0 



d d k 

(27T) 1 



7 e lkV exp 



-S4(Tc) 2 k a kf3M al3 



(71) 



where 



M a/3 = MX 



(72) 
(73) 



is a symmetric d x d matrix. Since we have the property 

V a M» = , 

we see that the matrix M a p has a zero eigenvalue corresponding to an eigenfunction in the T> direction: 

V a M aP = . (74) 
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In order to carry out the integral in Eq.(|71|), we must set up a coordinate system which singles out the T> direction 
and n orthogonal directions. Thus we construct an orthonormal set (V, t;^), for s = 1,2 ... ,n, which satisfy 



£i s) v a = o 



tWfV) — x , 



In Appendix C we show that we can write 



4 a) = E^ 

and it can be shown generally, see Appendix C, that 

(det A)(det A) 

where N is the n x n matrix 



1 



1 



det N V 2 



AT tVfV 

We then make the change of variables in the integral Eq. (|7l]) from k to 

n 

ka = ^2t v i^ +t d T> a 

v\ 

which clearly has a Jacobian of one. The integral of interest is then given by 

m- 



d d t j+ ,*n.i> 4+ £<y)ir ^ 4 fT \ 2 i + D 



(2ir 



(75) 



(76) 



(77) 



(78) 



(79) 



(80) 



(81) 



where v and v 1 range from 1 to n. One can then do the t^ integration to obtain a <5-function. The rest of the 
integrations are over gaussian fields governed by the matrix 

Quw = QH^Mce (82) 
which does not possess any zero eigenvalues. Evaluating the standard gaussian integral leads to the result: 



J(Q = 8(V-V)- 



-exp 



2Si{Tc) 



(27r5 4 (rc) 2 )"/ 2 {det Q)V2 
where the n-vector V v is defined by 

V — f^V 

Note, if we insert Eq.(]84|) into Eq. j83|) , we see that we need only the matrix 
and 



AO 



S(D ■ V) 1 



1 



; exp 



1 



2S 4 (Tc) 



;V a VfjR a i3 



(83) 



(84) 



(85) 



(86) 



(Tc) n (27r5 4 )™/ 2 {det Q) 1 / 2 ' 
We show in Appendix D that the matrix R, defined by Eq. (|S5|) , can be put into the very simple form 

R<xf3 = £,a£,/3 > ( 87 ) 

and finally J(£) is given by: 
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V ^ (Tc) n (27rS , 4 )™/ 21 1 
Inserting this result into Eq.(|68|) gives 

d n 



1 



25 4 (r c ) 



1 



(2ttS 2ttS 4 ) 



n/2 



n n^ 2 (o 



fl — l u—l 



5(V ■ V) 1 
(rc) n (2nS 2 ) nd / 2 ' 



where 



n d 



n d 



A ® = T EE^+aT^E E 



If we make the rescalings 



and 



then 



and we obtain 



n Q P[V] 



v 



•32 



P^(5 2 )" /2 P 



I I /5 2 / Sf W ' 



(27r) n (r c ) d V s 4 \s s 4 

where the dimensionless integral /(V) is defined by: 



I(V) 



d n 



m = I n n ^ ■ 

/2— 1 i/=l ^ ' 



with 



4>(o = E 
i/=i 

We see that the characteristic speed 



d d 

E(^) 2 + E 

a=l a, (3=1 



« = (rc) — 



(89) 



(90) 



(91) 



(92) 



(93) 



(94) 



(95) 



(96) 



(97) 



has been introduced into the problem. 

We can construct a form for the integral I(V) which does not involve a unit vector in the 5-function via the following 
rearrangements: 

5{V-V) = \V\S(V-V) 



1c 

YD 



5{V ■ V) 
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(detN) 



V 2 -, * 



(2tt)" /2 



V 2 6{V- V)e-* Z » N »■»'*•>' . (98) 



m = j n n rf c y 4 (o« • ^)^^ 



Inserting this result into the integral /[V] gives 

.1 " / ,]"- - , 

11 11 

/i — 1 I/— 1 

where 

VII. CASE n = 1 AND d = 2 

It is straightforward to work out P[V] and no for the simplest case of defect lines in two dimensions. The key 
simplifying aspect in this example is that the matrix reduces to a vector £ a and 

»»i(0 = WiU ( 101 ) 



with 



We have from Eq.(94) that 



V 2 = t 2 . (102) 



I I /5 2 / Sf \ 1/2 



nv) = Jd 2 ( f ^pf^K) ■ ^jtj e ^ A °' e) < 104 ) 

and 

M£,*) = & + V a Vf,£ a £fi + * > & ■ (105) 

Then, since the integral is isotropic, we can pick V to be in the x-direction and use 2?(£) • V = V£ y in the (5-function 
to do the integral over £ y and obtain 

_L_ fdi [ dz ^ ,-iCi^ 2 +v 2 ] 



^ = 7^ / d ^ / r-a75« e-^-li+* 1 . (106) 



The remaining integrals are elementary and we obtain the results 



and 



I(V) = —~ \ (107) 

(ttV) (1 + V 2 ) 2 



(2tt) Tc V S 4 V^o^4/ (ttV) (1 + l^ 2 ) 2 
We can then obtain the string density no in two ways. We can compute it directly from 
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n =< \V{$)\5{if) > 



(109) 



n<%;\v(o\G(0 (no) 



where G(£) is the integral over b of Eq.(|32|) given by 

= (2ttS' )«/2 (2irS 2 ) nd / 2€XP 
Restricting the analysis to n = 1 and d = 2 gives 



- 1 V(n 2 

25(2) Z^^' 



(111) 



"" > d ^vm^k e ~ 2ke • (112) 



These integrations are elementary with the final result: 

-°=Kif ■ <™> 

The second method for determining no, which serves as a check on intermediate steps in the calculation, is to integrate 
no-P[V] over all V and use the fact that P[V] must be normalized. The exercise is straightforward and leads to the 
same result for no. We then have the final result for the probability distribution 

2 1 1 

P[V] = ^ Vd{l + (V/v) 2 ) 2 (114) 
and it is easy to see that this agrees with Eq.(Q) for n = 1 and d = 2. 

VIII. CASE n = 2 AND d = 3 

We turn to the physically important case of strings in three spatial dimensions. In working out this case, the key 
observation is that 

1 



^a{0 — 2 eQ >MW2 e ^i,f2£/jl£p2 ' (us) 

is a quadratic form in £. This suggests that we use the integral representation of the 5-function 

8(P@ ■ V) = [ ? e ik ^ (116) 
J 2tt 

to write the integral of interest, Eq. (p9|) , in the form: 



3 2 



W'/nn^/ if/^(o.-««^ . .117) 

J fj,= l u=l v ' J J 

The matrix appearing in the gaussian is given by 

where we have introduced k a — kV a and have used the result 

-^(-^W)^=*'%) • (H9) 

We see that in principle we can carry out the £ integration since it involves a product of polynomials times a gaussian 
weight. Let us define the integral 



12 



mz) = / n n ■^^>~ kq ^' ( i2 °) 

and 

s 1 **) ■ <121) 

A significant simplification occurs if we realize that gradients with respect to k pull down factors of V and we can 
write 

L(k,z) = ViL (k,z) (122) 

where 

Lo(k,z)= / II II "Si' - '' •^^ • (123) 
/i— i i/— i ^ ' 

Thus we are left with a set of gaussian integrals. 

The gaussian integral can be carried out if we think of Q^, as a 6 x 6 symmetric matrix. If we can find the 
eigenvalues Aj, i = 1, 2 . . . , 6 of this matrix, then we have 

L (k,z) = ^L= . (124) 

VnliA* 

In Appendix E we discuss the relevant eigenvalue problem. The result of the analysis there is that the six-dimensional 
problem factorizes into a product of two three-dimensional problems. These three-dimensional eigenvalue problems 
reduce to cubic equations and the product of the three associated eigenvalues can be read off from the associated 
characteristic equation with the final result: 

£ (fc,z) = — = L =— 1 ^_ . (125) 

V (1 + z 2 ){l + z 2 + V 2 + k 2 ) + (k ■ V) 2 V (1 + z 2 ){\ + V 2 ) + k 2 + (k- V) 2 

We must then apply to L (k 1 z) and set k — kV. After a great deal of algebra one finds a complicated result for 
I[V] which still requires integration over k and z. It turns out that it is wise to first do the k integration. All of the 
contributions are proportional to integrals of the form 



/ 



2tt (1 + z 2 + V 2 k 2 )P+ 3 2ttV{1 + z 2 ) 5 / 2 V ' 

where p takes the values 0, 1 and 2 with k = K\ — jg, kq — g^fgy- The final integrals over z can all be expressed 
in terms of the integrals 



j - f l± I I M27) 

SUS2 J 27T ( l +z 2 + T >2 )si /2 ( l + z 2 )S2 /2 ^ Z <> 

which, for integer Si and S2, can be worked out analytically. After an enormous amount of additional algebra we 
obtain the very simple result 

I[V] = i = . (128) 

V (1+ T/ 2 ) 5 /2 



This leads back to the result 
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Integrating over all V we easily obtain the density of strings 

1 £ 2 
no = —w- • 

Putting this result back into Eq.( |129| ) gives the final result 

3 1 1 



P[V] 



which agrees with Eq.(pl) for n — 2 and d = 3. 



47T v 3 V (1 + V" 2 ) 5 / 2 



(130) 



(131) 



IX. WALLS IN THREE DIMENSIONS 

Let us conclude with the example of walls in three dimensions (n — 1, d — 3). We have from Section V that 

V 2 = (VV>) 2 (132) 



and the velocity of the wall is given by 



(133) 



Following the same path as for the point and string defects we find that the velocity probability distribution is given 
by 



n P[V] = J db d 3 mS (V + Tcgb/e) G(l b) 
where G(£, b) is given by Eq. (|62|) with n = 1 and d = 3 and 

n P[V] = [db d 3 ^\S (V + Tcgb/Z 2 



(134) 



1 



(2 7 rS*o) 1 / 2 (2^5 4 )!/ 2 (27r5 2 ) d / 2 
Again inserting the integral representation for the (^-function and doing the b integration, we obtain 



where 



n P[V] = 



1 



y/2nS (2nS 2 ) d 



d 3 m\exp 



— ? 
2S 2 



(2^) 



-l5 4 (r c ) 2 (fc.o 2 /e 2 



(135) 



(136) 



(137) 



The integral in Eq.(137) can be carried out if we introduce the orthogonal set of coordinates £, b^\ b^ and we choose 

i = 6 (1) x S« . (138) 

Then we obtain 



J(0 = 5{V -b w )8(V -6 (2) ) 



V27rS , 4 (rc) : 



25 4 (r c ) 



-(v-0 2 



(139) 



If we make the same change of variables given by Eqs.(BlJ) and (p2j) we find the result 



n P[V] 



1 1 



v 3 (2tt) 5 / 2 V & 



-J[f] 



(140) 
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I[V}= I d 3 Z?5tf-bW)6(f-bW)e H ?+(S ' )2 J . (141) 



and 



If we now define 

&W = x I* (142) 

and 

sw=e><s (2) (143) 

and assume that is in the z-direction we obtain 

I\V] = J d 3 tt 4 5(VbP)5(VbW) e -i{?+^ 2 ] . (144) 

The key point then is that 

5{VbP)8(VbW) = l-m,m y ) (145) 
v 1 

and the integral reduces to an elementary one-dimensional integral 

AV\ = i- [ dUte-^v*] _ (146) 



V 



Inserting the result of doing this integral back into Eq.(140) leads immediately to the result that P[V] is again given 
byEq.© . 

X. DETERMINATION OF v 

The determination of So, §2, S4 and v requires a theory for the auxiliary field correlation function 

C (12) = i(m(l)-m(2)) . (147) 
n 

If we follow the development in Refs.[ ^Jand [ we can write in the scaling regime, 

C (12) = S e- f2/{2L ^ (148) 



where f = f{ — r^. We do not need So, 5*2, and S4 individually. We only need the combination given by Eq.(|97|). We 
easily find 

o 2d , v0 Ted 



XI. CONCLUSIONS AND QUESTIONS 

The results of this work are simply summarized by Eq. (Q) . There remain several open questions. The most important 
question is whether Eq.(Q) for P\V] corresponds to the results found in the real world. It is highly desirable to measure 
P[V] for the rather broad range of systems covered by Eq. ([!]). 

How is this result for the defect velocity probability distribution changed at higher order in perturbation theory? 
So far, P{W) has only been calculated for the lowest order in the perturbation theory developed in Ref.( ||). Scaling 
arguments^ would indicate that the large velocity tails will not be modified at higher orders in perturbation theory. 
This question will be addressed soon. 

The results for P\V] given by Eq.(|l|) have only been proven for a subset of the range n < d. What about for 
n — d — 2 = 2? This question is somewhat academic since it is outside the range of physically accessible spatial 
dimensions. 
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APPENDIX A: V 2 

Clearly one of the important quantities entering into the discussion of the string velocity probability distribution 
is the quantity T> defined by Eq.(|l3|) and its square: 

This quantity cn be written in a simplified form if we realize that 

where / is a set of products of <5-functions giving all matched pairs between the unprimed and primed sets. There 
are minus signs if the matched pairs are an odd number of permutations of the labels in order to return to the order 
1,2, ... ,n. We have then 



det(N) . (A3) 



APPENDIX B: MATRIX M£ 

The matrix M£ defined by Eq.@ 

can be put into another useful form by inserting the explicit form for T> given 



by Eq.(pO|) and using the result given by Eq.(|A2|). An intermediate step gives 



T) — — P v ' x f- v '' 2 t"' n 



= -/[//, fx 2 . . . (j, n ; // x /4 . . . ^—fsy a ..^JXQ t • • • 

— -Zv' x v> 2 ...v' n ivi£,, ■ ■ ■ (Bl) 



Inserting this into the definition of gives 



1 1 



M v = — - - e P V2 £ Vn e * i ,f Vl t v ? £"V 



(B2) 
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This expression is particularly useful if we consider the quantity 

1 1 



Q M ^ — ^ 2 (n - iy ev ^---v n (i v' 1 u' 2 ...v' ri N u , u ' i N l , 2l / 2 . . . Af„„< 



7 ^ 1/2 ■ ■ - Vn ^-l/' ,1/2 ■ ■ -l^n 



V 2 {n - 1)! 

= -5u,v> ■ (B3) 

APPENDIX C: ORTHOGONAL COORDINATE SYSTEM 

We will construct an orthonormal coordinate system with the basis vectors 

n 

& 4) =X>„£ ■ (Ci) 
i/=i 

These basis vectors are orthogonal to V since 

%P a =Q . (C2) 

Since we require 

£i s) £k s,) =s ss > , (cs) 



we have immediately, on inserting Eq.(Cl) into Eq.(|C3j), that 

S S3 > = A SV A S , V ,N VV , . (C4) 

If we take the determinant we obtain 

1 = det A det N det A (C5) 



or, using Eq.([A3|), 

dpt. A dpt. A = 

det N £>2 



det A det A = — — = -4- . (C6) 



We can obtain a realization of the matrix A sv by constructing the £a directly using the basis set 



X 



(s) = e 2— /"^ ( C7 ) 



which are not orthonormal. However they can be constructed to be orthonormal using the Gram-Schmidt orthogo- 
nalization process. Thus each is a linear combination of the Xa which is proportional to a linear combination of 
Thus one can extract an explicit expression for the matrix A. This explicit expression is not needed here. 



APPENDIX D: INVERSE OF MATRIX Q 

We need the inverse of the matrix 



If we use the expression given by Eq.(|Cl|) for the basis vectors we obtain 



Qu„> = ^ig'^^ . (Dl) 
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Qw> = A vV HA v , v ,^ Ml Ml . (D2) 



We then use the identity Eq.(B3) twice to obtain 
This immediately tells us that 

det Q = det A det A = -tr- . (D4) 
V 2 

We still need to construct the inverse of Q. It is easy to see that this is given by the product 

Q- 1 = A' 1 A- 1 . (D5) 

However we only need the elements 

R a0 = & (Q~XJP 



m ■ (D6) 



APPENDIX E: EIGENVALUE PROBLEM 



We need to find the six eigenvalues for the matrix Q defined by Eq.(118). The analysis can be carried out by looking 
at the action of Q when acting on the six basis vectors: ipi — V^z Vl ip 2 — V^c v , ips — b^z Vl ip^ — b^c u , ip^ — k^z V) ip^ — 
k^Cv where we have introduced c v — e UV 'Z u i, and b — k x V. The action of Q acting on these states is given by 



Qipi = {1 + z 2 + V 2 )^ 1 +ikip i 


(El) 


Q^ 2 = (1 + V 2 )^ 2 - ik^ 3 


(E2) 


Qip 3 = (1 + z 2 )ip 3 ~ iki> 2 + ik ■ Vtpe 


(E3) 


Qip4 = ip4 + iktpi — ik ■ y-05 


(E4) 


QVs = (1 + z 2 )t/j 5 + V 2 k ■ ViPi 


(E5) 


Q^j e = + V 2 k ■ V4> 2 ■ 


(E6) 



We see that these equations decouple into two cubic systems (ipi, ip4, ips), and (4>2, V"3i "06)- We then only need the 
products A1A2A3 and A4, A5A6. It is well known (see Ref.[ [35|) that if one has a cubic characteristic equation 

A 3 + a 2 X 2 + aiA + a = , (E7) 

then the product of the three roots is given by 

A1A2A3 = -ao . (E8) 
In the present case we easily find for the two sets: 

a ( o } = -(1 + z 2 )(l + z 2 V 2 + k 2 ) - (k ■ V) 2 (E9) 

and 

a (2) = -{l + z 2 )(l + V 2 )ls-+k 2 )-{k-V) 2 . (E10) 
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